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Newly born magnetars are good candidate sources of ultrahigh energy cosmic rays. These objects 
can in principle easily accelerate particles to the highest energies required to satisfy the ultrahigh 
energy cosmic ray scenario [E ~ ]^q20-2i gy)^ thanks to their important rotational and magnetic 
energy reservoirs. Their acceleration mechanism, based on unipolar induction, predicts however a 
hard particle injection that does not fit the observed ultrahigh energy cosmic ray spectrum. Here 
we show that an adequate distribution of initial voltages among magnetar winds can be found to 
soften the spectrum. We discuss the effect of these distributions for the stochastic gravitational wave 
background signature produced by magnetars. The magnetar population characteristics needed to 
fit the ultrahigh energy cosmic ray spectrum could lead in most optimistic cases to gravitational 
wave background signals enhanced of up to four orders of magnitudes in the range of frequency 
1 — 100 Hz, compared to the standard predictions. These signals could reach the sensitivities of 
future detectors such as DECIGO or BBO. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The origin of ultrahigh energy cosmic rays (UHECRs) 
is still unknown (see [l| for a review). From the point 
of view of simple particle confinement energetics (the so- 
called "Hillas criterion" Q), young millisecond magne- 
tars are one of the most promising candidate sources. 
Magnetars are isolated neutron stars with extremely 
strong surface dipole fields of order ~ 10^^ G and fast 
rotation at birth with initial rotation period Pj ~ 1 ms 
(see for reviews). Because the source of energy for 
their radiative emission is magnetism, their dissipative 
properties are distinct from those of radio pulsars. Their 
existence was postulated by Ref. and they are ac- 
cepted as a plausible explanation for Soft Gamma Re- 
peaters and Anomalous X-ray Pulsars (e.g., [THSj]). 

Ordinary pulsars have long been discussed as good can- 
didates to accelerate charged ions through unipolar in- 
duction mechanism (see, e.g., (loj and references therein). 
Rapidly rotating neutron stars generally create relativis- 
tic outflows ("winds"), where the combination of the ro- 
tational energy and the strong magnetic field induces an 
electric field E = — v x B/c (where v and B are the ve- 
locity and the magnetic field of the outflowing plasma). 
The wind thus presents voltage drops where charged par- 
ticles can be accelerated to high energy. Ordinary pulsars 
however do not supply enough energy to reach the high- 
est energies {E > 10^° eV). 

Magnetars on the other hand possess important rota- 
tional and magnetic energy reservoirs at birth that should 
enable them to accelerate easily, in principle, particles to 
E > 10^° eV [m, They were introduced as possi- 
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ble progenitors of ultrahigh energy cosmic rays (UHE- 
CRs) during the "AGASA era" . If cosmic-rays originate 
from cosmological distances, their flux at the highest en- 
ergies {E > 6 X 10^^ eV) should be suppressed due to 
interactions with the cosmological background photons, 
crea ting a feature in the spectrum called the "GZK cut- 
off" [11 [H]. The fact that the AGASA experiment did 
not observe such a feature led Ref. [ll| to develop a 
model of acceleration of ultrahigh energy iron nuclei in 
youn g st rongly magnetized Galactic neutron star winds. 
Ref. hj followed the same trend, but suggested that 
the hard injection spectrum produced by each magne- 
tar could account for the absence of GZK cut-off, even 
with an extragalactic magnetar population scenario. 

The latest experiments report however that a suppres- 
sion reminiscent of the GZK cut-off is present at the high- 
est energy end of the UHECR spectrum [H, The 
cosmic ray spectrum observed by the Pierre Auger Ob- 
servatory can be described as a broken power- law, 
with spectral index A 3.3 below the break (called "an- 
kle") around 10^^'^ eV, and A ~ 2.6 above, followed by a 
flux suppression above ^ 10^^'^ eV [l3|- The mainstream 
UHECR models view the ankle region as a transition be- 
tween cosmic rays produced by Galactic and extragalac- 
tic source populations. Given this situation, the hard 
injection spectral index of A = 1 produced in the above 
magnetar scenarios is no longer an advantage. Figure [T] 
shows indeed that it is challenging to fit the observed 
spectrum down to the ankle energy with extragalactic 
sources injecting such a hard spectrum. The introduc- 
tion of various source emissivity evolutions (described in 
Section IHI[) are insufficient to reconcile the calculated 
and observed spectra, even for strong evolution cases. 

In this paper, we show that an adequate distribution of 
initial voltages among extragalactic magnetar winds can 
be found to soften the overall UHECR spectrum. These 
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distributions result in some cases in tighter constraints on 
the magnetar population rate that is required to account 
for the observed UHECR flux, than suggested in Ref. [13] . 

The magnetar model envisages that a wound-up, 
mainly toroidal magnetic field with strength Bt > 10^^ G 
characterizes the neutron star interior |18| . Such strong 
internal fields should lead to a substantial deformation 
of the neutron star and thus to the emission of gravita- 
tional waves, provided that the magnetic distortion axis 
and the rotation axis of the star are not aligned 0, [l^ . 
Young millisecond magnetars should thus be strong grav- 
itational wave emitters [2ll - [23l |. 

Gravitational wave signals from individual magnetars 
might be detected by instruments of the generation of 
Advanced LIGO [1^, and reveal properties about their 
magnetic field and initial rotation period that are crucial 
to probe magnetars as UHECR accelerators. However, 
such signals should not be observed in coincidence with 
UHECR events. Indeed, the time delay that cosmic rays 
experience by magnetic deflection during their propaga- 
tion in the intergalactic medium, relative to the gravi- 
tational waves going in geodesies, is more than several 
hundreds of years. 

Another gravitational wave signature that one can seek 
is the stochastic background emitted by the ensemble of 
magnetars. The signal calculated up to now assuming a 
population of magnetars with identical properties, were 
weak especially in the frequency range that should be ob- 
served by future experiments such as BBO or DECIGO 
[23 |. We demonstrate in this paper that the distribu- 
tions of magnetar characteristics required to account for 
the observed UHECR spectrum could, in some cases, sig- 
nificantly enhance the signal level in the frequency range 
1 — 100 Hz. If such signatures were detected, they could 
help determine if magnetars are indeed capable of accel- 
erating the highest energy particles in the Universe. 

The layout of this paper is as follows. In Section|TTl we 
provide a synthetic picture of the production of UHECRs 
by newly-born magnetars, not entering into the detailed 
modeling of the acceleration mechanism or the escape of 
particles. In Section IIIIl we calculate the UHECR flux 
obtained for various magnetar distributions and discuss 
the magnetar occurrence rate necessary to account for 
the observed cosmic ray flux. We calculate the implica- 
tion of such distributions on the stochastic gravitational 
wave background in Section IIVI Our results are further 
discussed in Section Ivl 



to be further investigated on several issues, for instance 
on the nature of the ions injected in the wind, on the 
mechanism through which the current in the wind taps 
the available voltage (i.e., the actual acceleration mecha- 
nism) , and on the escape of the accelerated particles from 
the wind, and the surrounding supernova envelope psl j. 
Discussions on these subjects can be found in |12l |. 



A. Maximum acceleration energy 

The internal, mainly toroidal magnetic field par- 
tially threads the magnetar crust, making up a mainly 
poloidal magnetosphere with surface dipole strengths 
Bd ~ 10^^ G that are req uired to account for the ob- 
served spin-down rates llq. This dipole component de- 
creases as B{r) — {l/2)Bd{R^/r)^ according to the dis- 
tance from the star's surface r, with the radius of the 
star. Beyond the light cylinder radius i?L = c/fl, the 
dipole field structure cannot be causally maintained and 
the field becomes mostly azimuthal, with field lines spi- 
raling outwards and with strength decreasing as B{r) ^ 
B{Ri^){Ri^/r). The out-fiowing relativistic plasma at 
r > _Rl (the magnetar "wind") thus has magnetospheric 
voltage drops across the magnetic field of magnitude 

mm 

$wind - rB{r) ^ i?LB(i?L) - ^ (1) 

= 3 X 102^33^^4 V , (2) 

where Q.^ = 17/10^ s~^ is the angular velocity of the 
star and /i33 = /i/lO'^'^cgs its dipole moment with fi = 
BdRl/2 = 10^3 cgs(Bd/2 X lO^^ G)(i?*/10 km)^. 

Assuming that particles with charge q experience a 
fraction r? o f this voltage drop, they should gain the en- 
ergy ^H^: 

Bin) = gry$wind - - 3 X lO^iZryi 17^^33 eV , (3) 

where we define rji = rj/OA. For magnetars beginning 
their lives with millisecond rotation periods, particles in 
principle can achieve ultrahigh energies of E > 10^° eV 
by acceleration in the wind. 



B. Injection spectrum 



II. MAGNETAR ENERGETICS FOR UHECR 
ACCELERATION 

In this section, we provide a synthetic picture of the 
production of UHECRs by newly-born magnetars, fol- 
lowing the work of Refs. [HI, These authors propose 
that UHECRs are accelerated in the relativistic wind of 
rapidly spinning and strongly magnetized magnetars by 
unipolar induction. The toy model described below needs 



Inside the light cylinder, the magnetosphere corotates 
with the star and the ion density corresponds to the max- 
imum current density of particles extracted from the star 
surface, i.e., the Goldreich- Julian charge density pcj ^26, ]. 
Assuming that this current is entirely tapped in the wind 
for acceleration, one can write the instantaneous particle 
injection rate N-^ as a function of and thus of the parti- 
cle energy at a given time E. The energy spectrum of the 
particles accelerated by a magnetar during its spin-down 
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reads: 



dNj 



dfl dE 



(4) 



The spin-down of a magnetar is driven by electromag- 
netic energy losses and gravitational wave losses fsee for 
example Section 10.5 of Shapiro & Teukolsky [lO])- Ex- 
pressing these losses in terms of ft and thus E, Eq. ^ 
can be transformed into 1131: 



dA^i 

Je 



9 c'l 



AZe^iE 



E 



(5) 



where the critical gravitational energy at which gravity 
wave and electromagnetic losses are equal reads: 



E, 



5 Zrjeii^ 



3 X 10 



20 



eV. 



(6) 



with = //lO'^^gcm^, the principal moment of iner- 
tia of the star. The ellipticity of the magnetar, created 
by the anisotropic pressure from the interior magnetic 
field, can be evaluated numerically by: e ~ 10^^ — 
{Bf))/{4 X 10^^ G), where the brackets denote a volume 
average over the entire core (see, e.g., [l9l \2^). 

We note 62 = e/10^^. We neglect the influence of r- 
mode instabilities on the magnetar spin-down: this effect 
should mainly affect the cut-off of the injected energy 
spectrum by modifying the energy loss component due 
to gravitational wave losses ^12]. 

Equation ([S]) gives the energy spectrum of cosmic rays 
that are injected by a single magnetar in the interstellar 
medium, over 1 — 2 hours to go down to ankle energies 
(see Eq. 28 and 29 of Ref. [l3| for evaluations of the spin- 
down time). 



III. INTEGRATED UHECR FLUX FROM 
MAGNETAR POPULATIONS 

A. Uniform distribution 

We calculate the flux of UHECR on Earth produced by 
a population of magnetars. We first assume that all niag- 
netars have identical physical properties and that they 
occur at a rate rim per unit volume per year in our lo- 
cal Universe. At energies above the ankle, the influence 
of the intergalactic magnetic fields on the shape of the 
UHECR spectrum should be negligible. The cosmic ray 
spectrum dN/dE then follows the relation: 



d r ■ dN 

I ^ 

dE V dE 



dN, 

= Wgoom "m -777 ■ (7) 



Here VFgcom is a geometrical factor that accounts for the 
fact that all magnetars cannot inject ions from the stars' 
atmospheres into the wind in the rotational equator (see 
section 4.4 of Ref. For numerical applications, we 
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FIG. 1: Cosmic ray spectra for an uniform distribution of 
initial voltages among magnetar parameters (Eq. |9|. A pure 
proton composition is injected with the parameters chosen in 
the numerical apphcation of Eq. ([6]), with initial rotation ve- 
locity ^1,4 (i.e., Bmax = 3 X 10^^ eV and Eg = 3 x 10^° eV). 
The flux calculated analytically for no source evolution his- 
tory (black solid) is compared to the numerical propagation 
computations with difi'erent source evolution models. See text 
for the magnetar rate rim required to flt the data in each case. 
The spectra observed by HiRes and the Pierre Auger Ob- 
servatory are overlaid. 



will choose in the rest of this paper the value Wgcom = 
0.5. Integrating equation ([7]) using E = —E/T\oss (en- 
ergy losses due cosmological expansion and to interac- 
tions with the cosmological photon backgrounds) and 
Eq. ©, we get [i2: 



J{E) 



c dN 
4^'dE 



9 Jc3 



In 



£'max 1 + {E/Eg) 



E l + (S,nax/Sg) 



(8) 



(9) 



Emux is the maximum acceleration energy corresponding 
to the initial rotation velocity: i^max = E{ni). We as- 
sumed in this calculation that the magnetar birth rate 
njjj remains constant throughout time. Though mag- 
netars are bursting UHECR sources, the spread in their 
arrival time on Earth induced by deflections on the in- 
tergalactic magnetic fields should be sufficient to account 
for the continuous detection of particles (see discussion 
by Ref. [27] for the case of gamma-ray bursts, and Sec- 
tion IIV Al of this paper for estimates of the time delays 
due to deflections). 

Figure [1] presents the UHECR flux calculated analyti- 
cally following Eq. for magnetars with the parameters 
chosen in the numerical application in Eq. ([5]), with ini- 



4 



tial angular velocity fli^4 (i.e., i^max = 3 x 10^-^ eV), and 
injecting a pure proton composition (black solid line). 
The calculation of the term T\ciss{E) due to the cosmo- 
logical expansion and to energy losses on the cosmolog- 
ical photon backgrounds takes into account photo-pion 
production and pair production processes on the cosmo- 
logical microwave background (CMB) and the infrared, 
optical and ultraviolet background photons modeled by 
Ref. dl]. The rate of magnetars required to fit the ob- 
served spectrum is nm,o = 7 x 10~* Mpc""^ yr~^, which 
corresponds to 3.5% of the whole magnetar population, 
if one assumes that a magnetar birth rate of 10"'' per 
year per galaxy, and that the average galaxy density is 
of 2 X 10~^ Mpc"'^ 29]. One can note however that both 
HiRes [30] and the Pierre Auger Observatory [l3l report 
systematic uncertainties of order 20% on the absolute en- 
ergy scale of the spectrum, which should be considered 
for the evaluation of n,n,o- 

We also present the spectra calculated numerically us- 
ing the cosmic ray propagation code developed in [sTI - 
ISSj . assuming that each source injects particles follow- 
ing Eq. ([5]), with the same magnetar physical parame- 
ters as in the analytical case. Cosmological expansion 
losses, as well as losses due to interactions of UHECRs 
with the CMB and the infrared, optical and ultraviolet 
background photons modeled by Ref. [1^ are included. 
We examine the effects of three typical source emissivity 
evolution cases: (i) no evolution, (ii) the source evolu- 
tion follows the star formation rate normalized to unity 
at z = derived in Ref. |34] (labeled SFR): 



Rsfr{z) = 



1 + 7Mz 



1 + (z/3.3) 



5.3 



(10) 



and (iii) a strong source evolution case that for exam- 



ple Faranoff- Riley type II galaxies might follow 35 1 

36| 



These evolutions are described in detail in Ref. 
The uniform case is, as expected, in good agreement 
with the analytical calculation. The spectra are fit- 
ted by eye to the data and the magnetar birth rate at 
z — required in each case are: rim.SFR ~ 0.8 rim, o ~ 
5.6 X 10~* Mpc~^ yr"^ for the star formation type evo- 
lution, and nm, strong ~ 0.7 rim, ~ 5 X 10~^ Mpc~'^ yr^^ 
for the strong evolution. 

The spectra calculated for an uniform distribution of 
initial voltage drops among magnetars only mildly fits 
the observed UHECR spectrum. One would need another 
type of source to cover the energy range ioi^~i9-5 eV, a 
scenario which is not very attractive in terms of simplicity 
and fine tuning. We propose in the next sections to in- 
troduce a distribution of the initial voltage drops among 
magnetars to reconcile the observed spectrum with the 
magnetar scenario. 



B. Distribution of fii 

Equation © assumes that all magnetars have the same 
initial voltage $i = il?/i/c^. It is likely however that this 



voltage differs from star to star, because of different ini- 
tial angular velocities ilj, and/or differences in the mag- 
nitude of the surface dipole fields and hence on the dipole 
moment /i. The distribution of both these quantities 
among the magnetar population is basically unknown. 
One might then relax the previous assumption and in- 
troduce a distribution of magnetar birth rates according 
to the starting voltage, assuming a power-law: 



1 



drir, 



d$i «'i,max (*i,max/4'i,min)' 



1 



(11) 

with $i,min < 'fi < $i,max- As a function of the initial 
acceleration energy Ei, we get: 



dnin drim d'l'i 
where we defined: 

1 

X 



Ei 



1 



{E 

i,max/ ^i,min)^ ^ 1 

Equation ([9)) is then transformed into: 



J{E) 



dEi 



■dEi . 



(12) 



(13) 



(14) 



This integral leads to different results according to which 
parameter fij, fi, or both is distributed among the mag- 
netar population. For a distribution of initial angular 
velocities ili with a fixed dipole moment fi among the 
stars, the integral is calculated in Appendix El One can 
note that the slope of the overall UHECR spectrum be- 
low the high energy cut-off happens to be s, the slope of 
the distribution of ilj among magnetars. 

The resulting cosmic ray spectra when distributions of 
fli are included, for a pure proton injection, i?i.min = 
3 X IQis eV, £:i,max = 3 X lO^i eV and E^^ = 30, 300 
and 3000 EeV are presented in Fig. [2j The analytical 
results from Appendix [X] are plotted in solid lines. The 
results of our numerical particle propagation are shown 
in dotted lines, with colors corresponding to the critical 
gravitational energy Eg labeled for the analytical cases. 
For each Eg, the spectrum obtained for no evolution, SFR 
and strong evolution models are shown (from bottom to 
top in the low energy end) . The numerical calculation is 
performed by assuming that all magnetars inject particles 
according to Eq. (O, up to a maximum energy that is 
sampled according to Eq. (fT2|) . The agreement between 
the simulated and the analytical spectra in the absence 
of source evolution is good. 

Note that in principle, the value of Eg could also de- 
pend on the angular velocity f2i through the ellipticity e, 
as we will see in the next section. This dependence is not 
trivial as we will discuss, and we assume for simplicity in 
this section that Eg is fixed for all magnetars. 

Table [J recaps the values of the magnetar population 
parameters required to fit the observed UHECR spec- 
trum for the various scenarios represented in Fig. [21 
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FIG. 2: Cosmic ray spectra when distributions of fii are in- 
cluded, for a pure proton injection, i5i,min = 3 x 10^* eV, 
.max — 'J xlO^^ eV, M33, I45 and 771. Colors indicate Eg = 30, 
300 and 3000 EeV respectively. Analytical results from Ap- 
pendix |A] are plotted in solid lines. The results of our numer- 
ical particle propagation are shown in dotted lines, for no, 
SFR and strong evolution models (from bottom to top in the 
low energy end). The magnetar population parameters used 
to best fit the observed data are presented in Table U 



FIG. 3: Cosmic ray spectra calculated analytically when dis- 
tributions of Sli, fi and /i = f{Q,i) are included (different 
colors), for a pure proton injection, £"1,111111 = 3 x 10^* eV, 
iSi.max = 3 X 10^^ cV, fi33, I45 and rji. For the distribution of 
r^i, cases for Eg — 30, 300 and 3000 EeV are represented with 
decreasing thickness. The magnetar population parameters 
used to best fit the observed data are presented in Table [III 



The required magnetar birthrate is fairly high, rang- 
ing from 10 — 30% of the total supernova rate (of or- 
der IQ-'' Mpc-3 yr-i). The rate of magnetars is widely 
uncertain, but such values are advocated by some au- 
thors [131 • One may note furthermore that these rates 
includes mostly neutron stars that are are born with rel- 
atively low rotation periods (fii^min ~ 300 s~^, for /X33). 
The required rate is slightly lower when a source evo- 
lution is included. The values indicated in Table U are 
approximative, as the fits are done by hand and because 
of the intrinsic uncertainties in the observed data. Vary- 
ing i?i.max has little impact on rtm.o and s because most 
of the population contributes at low energies due to the 
decreasing power-law distribution. 

Finally, Figure [5] demonstrates that the three values 
of Eg represented are allowed by the observed data. 
Low and strong gravitational wave emission models 
are equally possible at the present experimental stage. 
Stronger statistics in the highest energy end of the spec- 
trum would help constrain the shape of the cut-off. 



weaker /j, will thus lead naturally to stronger fluxes and 
varying over [/J.i,min, Mi.max] softens the integrated cos- 
mic ray spectrum without requiring higher source densi- 
ties to account for the low energy flux. 

The critical gravitational energy Eg also depends on 
fi (Eq. [6|). However, this dependency is not straightfor- 
ward as the ellipticity e is also expected to depend on the 
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C. Distribution of fj, 

We now assume that fij is fixed and that the dipole mo- 
ment fi has a distribution among magnetars that lead to 
the distribution in voltages described by Eq. pT|) . This 
case is interesting, as the injection spectrum of single 
magnetars themselves scales inversely to /i (Eq. [S]). A 



TABLE I: Parameters for distribution of Q,i {Sec- 
tion IIII B)) . Magnetar density at z — 0, Un-i.o, and spectral 
indices s of the distribution of initial voltage drops required 
to fit the observed UHECR spectrum are indicated for vari- 
ous source model evolutions and critical gravitational energy 
Eg. We assume a pure proton injection, i5i,min = 3 x 10^* eV, 
£i,max = 3 x 10^^ eV, fj,33, I45, and r?i. 



6 



magnetic field. For magnetic field induced deformation 
of the star, assuming that the star interior is a compress- 
ible perfect fluid, one can consider that the ellipticity can 
be written (iq} : 

e = /3^A^^ . (15) 

Here, /3 is the magnetic distortion factor introduced by 
Ref. [19], which measures the efficiency of the interior 
magnetic field in distorting the star. This factor depends 
on the equation of state of the star interior and on its 
magnetic field geometry. Ref. [l^ find that the value of 
/3 can range between 1 — 10 for perfectly conducting inte- 
riors (normal matter), 10 — 100 for type I superconductors 
and can reach > 100 for type II superconductors. The re- 
lation above might not stand however when the interior 
magnetic field is dominated by the toroidal component 
[23 |. and the magnetar internal field being already close 
to or exceeding the critical field value (of order 10^^ G 
[sst ) , it is possible that no significant enhancement of the 
magnetar deformation happens beyond e ^ 10^* — 10^'^. 

In the following, we will rely on this relation for sim- 
plicity. The influence of Eg (hence of the dependency 
between e and fj,) is only noticeable at the very high en- 
ergy end of the UHECR spectrum {E > few 10^° eV), 
with variations that could not be measured with the ob- 
servational statistics of current detectors. Differences in 
Eg have an impact on the normalization that are be- 
low the systematic uncertainties of the observed UHECR 
spectrum. 

The analytical spectrum for a distribution of /i among 
magnetars is calculated in Appendix |B] In this case, the 
spectral index of the UHECR spectrum before the high- 
energy cut-off is given by s + 1. The results are presented 
in Fig.Ofor a pure proton injection, no source evolution, 
f^i = 10^ s-\ Si,„,in = 3 X 10i« eV, Si,„,ax = 3 X 10^1 eV 
and /3 — 700, which corresponds to e ^ 10^^. The corre- 
sponding values for Tim.o and s are presented in the third 
row of Table |TT1 The shape of the spectra is similar to 
that obtained for a distribution of ilj and superimpose 
well in the low energy region. This is explained by the 
fact that the spectra calculated in Appendix [X] and iBl 
have nearly the same expression, the main difference re- 
siding in the normalization. With lower values of /3, one 
recovers the shape of the spectra calculated for high Eg in 
the previous section. In Table HIl we have also indicated 
the normalization and spectral index values required to 
fit the observed spectrum for /3 = 200, 700 and 2000. 
These values correspond to Eg - 3000, 300 and 30 EeV 
respectively, for 145 and ^33. One can notice that the 
influence of /3 on rim^o f^nd s is very mild. 

Here, the magnetar birth rates (in Table lll| needed to 
fit the observed UHECR spectrum are lower of two orders 
of magnitude compared to the case treated in the previ- 
ous section. The rates are close to those calculated for 
the uniform distribution. This case is less tight in terms 
of population density as compared to the distribution in 
rotation velocities, but necessitates that all the magne- 



tars considered are sub-millisecond rotators at birth. It 
is also "disappointing" in terms of gravitational wave sig- 
natures, as it should lead to faint signals, as we will see 
in Section [D2l 
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TABLE II: Parameters for different distributions. Mag- 
netar density a.t z — 0, rim,o, and spectral indices s of the dis- 
tribution of initial voltage drops required to fit the observed 
UHECR spectrum in absence of source evolution are indicated 
for different distribution models, critical gravitational energy 
Eg, distortion parameter j3 and a as defined in Eq. (|16[) . 
We assume a pure proton injection, i5i,min = 3 x 10^* eV, 
.max — 3 X 10^^ eV, I45, and rji. 



D. Distribution of jj. — /(Jli) 

It is plausible that the surface dipole magnetic field 
-Bd and its initial angular velocity Jlj are not independent. 
Ref. Jg] conjectures that an efficient aw-dynamo operates 
during the formation of the magnetar. For proto-neutron 
stars born with initial periods Pi, one can calculate that 
the saturation magnetic field (when energy equipartition 
is reached between the fiuid and the field) on small scales, 
generated by kinematic growth of seed fields under dy- 
namo action, is of order Bgat ^ 3 x 10^*" G (1 ms/Pj). 
The observed large scale dipole field Pd can be regener- 
ated from this field differential rotation and convection 
processes. One can thus assume Pd to be a constant 
fraction of the saturation field and write the following 
relation [s^: 

= «— , (16) 

TT 

where we will choose arbitrarily a G [10^^, 10^'^] Gs. The 
maximum value of a is chosen to account for the order 
of magnitude of the strongest dipole fields suggested to 
exist in magnetars, assuming fti — lO'* s"^. 

In reality, it is likely that the relation above be mod- 
ified, for example if the ratio between Pd and Pgat is 
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a function of Pi (see the corresponding discussion in 
Ref. ^]). Because these dependencies are not eluci- 
dated, and for the sake of simphcity, we will stick for 
the time being to the above formula. 

The situation is then similar to the previous section: 
magnetars birthrates are distributed according to fi, with 
a different relation between (or Ei) and /i. The angular 
velocity as a function of reads: 



o 2^ 



The initial voltage can then be expressed as: 



^ ^ 3 

qr] 



with A 



qri 



(18) 



The analytical calculation of the spectrum is then 
straightforward and can be found in Appendix [C] We as- 
sume in this calculation that E^^ is fixed among magnetars 
to avoid a complicated cut-off function at the highest en- 
ergies that would not affect strikingly the normalization 
of the spectrum, nor the spectral index of the distribution 
of initial voltages as discussed in Section IIIICl Also be- 
cause, as discussed in the same section, the dependency 
of £ on /i is not clearly determined. 

As expected, the spectrum found matches perfectly the 
one calculated for a distribution of fii (see Fig. The 
value of nni.o scales as a~'^/^ . 

In Table we calculated the values of the normal- 
ization and the spectral index needed to fit the observed 
spectrum in absence of source evolution. The cases for 
SFR and strong source evolutions can be roughly ex- 
trapolated through the following relations, as we checked 
in Section IIII Al and with supplementary simulations: 

"-m,0,SFR ~ 0.8ri,„^o,noevol, and SSFR '~ Snocvol — 0.2, 

for the star formation type evolution; nm^o^strong ~ 
0.7nm,o,nocvoi ~ 5 X 10"* Mpc"^ yr"\ and Sstrong ~ 
Snocvoi — 0.4 for the strong evolution. 



IV. IMPLICATION FOR THE DIFFUSE 
GRAVITATIONAL WAVE SIGNAL 

A. Gravitational waves from a single magnetar 

Rotating neutron stars are expected to emit copious 
amounts of gravitational radiation, mainly at its rotation 
frequency and at twice its rotation frequency, provided 
that it deviates from axi-symmetry ( [401] for a review). In 
the case of a magnetar, the deviation from axi-symmetry 
is believed to be caused principally by the neutron star's 
internal magnetic field (19i] . The energy losses from grav- 
itational waves for a star with ellipticity e can be esti- 
mated as follows [41II: 



E, 



32 GPe^Vfi 



(19) 



The gravitational spectral energy emitted by a single 
source between frequency v = ri/(27r) and v + Av thus 
reads: 



dE, 



gw 



Av 



2n 



K 



(20) 



with iye[0, i^i] ,(21) 



(17) 

where i>i Qi/ir is twice the initial spin frequency and 



K EE 



2887r4 GPe^ 72tt^ /S'^Rt 2 



5 c^GI 



(22) 
(23) 



The last equality assumes that the deformation of the 
magnetar is due to its internal magnetic structure, and 
that it is a quadratic function of the amplitude of the 
magnetic dipole moment, as in Eq. (|15|) . 

The detectability of gravitational wave signals from 
individual magnetars has been investigated by several 
authors [H, |4§| who find that the signals could be de- 
tectable by Advanced LIGO-class detectors up to the 
distance to the Virgo Cluster. Though such a detec- 
tion would be a watershed in several fields of high-energy 
astrophysics, note that gravitational radiation from sin- 
gle magnetars cannot be observed in coincidence with 
UHECR events. The delay induced by extragalactic mag- 
netic fields of mean strength B and coherence length 
on particles of charge Z and energy E with respect to 
gravitational waves over a distance D reads [H, : 



6t ~ 2.3 X lOVrs^^ 



D 



lOMpc 
V 1020eV 



A 



B 



0.1 Mpc, 

2/^x2 



10-9 G 



(24) 



It is likely that magnetic fields at this level of ^ nG are 
present in our local supercluster, over a few megaparsecs. 
For homogeneous intergalactic magnetic fields of lower 
overall strength {B < 10"^^ G), the time delay could be 
shorter than a year over 100 Mpc. However, the crossing 
of one single magnetized filament (size fi , field strength 
B and coherence length Ai) will lead to a slight deflection 
that will induce a time delay with respect to a straight 
line of order 



5ti ~ lO^'yrsZ^ 



2 Mpc/ V 0-1 Mpc 



A,, 



E 



1020 eV 



B 



10-8 G 



X (25) 



(26) 



For the same reason, it is unlikely that even for magnetars 
inside our Galaxy, gravitational wave signals could be 
detected in coincidence with UHECRs. 
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B. Diffuse gravitational wave signals expected for 
magnetar populations that are sources of the 
observed UHECR flux 

The possibility of detecting gravitational wave back- 
grounds generated by an ensemble of neutron stars and 
of magnetars more specifically, has been discussed by sev- 
eral authors [47 - 51.] . The background results from the su- 
perposition of the signals of the population of emitters, 
integrated over the whole history of their evolution. 

The spectrum of the gravitational stochastic back- 
ground can be characterized by the quantity r2g„(i'o), de- 
fined as the present-day energy density per logarithmic 
frequency interval, in gravitational waves of frequency 
vq, divided by the critical energy density of the Universe 
pc^. For a population of sources, each emitting gravita- 

this quantity can be 



tional waves according to Eq. (|2T 
expressed [5^.l53|: 



10" 
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/ , LIGOIII/ET 




DECIGO Adv. 
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ngw(^o) = 5.7 X 10"^^ 



0.7 

ho 



1^0 

Rcvoi{z) dEg^ 



■ di^i X 



lin ^'^i 

1^0(1 + 2)] dz, (27) 



where i/q and v = vo{l + z) are the frequencies in the 
observer and in the source frame respectively, U{z) = 
[SlA-l-Om(H-z)'^]^/^, i?evoi(^) IS the dimensionless source 
evolution rate normalized to 1 at z = 0. The integral 
upper bound Zsup{i^i) is given by: 



if t/Q < — — 
— — 1 otherwise, 

I'D 



(28) 



where Zmax = 6. A distribution of the source density 
according to the initial frequency Vi (corresponding to a 
distribution in "I>i) has been introduced in the calculation 
of the background signal. 

The number of magnetars is large enough for the time 
interval between events to be small compared to the du- 
ration of a single event. One can show indeed that the 
duty cycle 



A(z) = y n„i^o-Rovoi(z')T-gw(l + 2')-^(^')d2' 



(29) 



defined as the ratio between the duration of the events 
and the time interval between successive events satisfies 
the condition A ^ 1. Tgw is the average duration of a 
signal produced by a magnetar; the assumption that this 
duration is of the order of the magnetar spin-down time 
leads to A > 10'^. Their gravitational wave signal can 
thus be considered as a continuous background. Such 
backgrounds obey the Gaussian statistic and are com- 
pletely determined by their spectral properties. 

The optimal strategy to detect these signals is to per- 
form a correlation between two or more detectors, possi- 
bly widely separated to minimize common noise sources 



FIG. 4: Energy density of the stochastic gravitational wave 
background f2g„ produced by magnetars as a function of the 
observed frequency, assuming a distribution of fii among mag- 
netars, as in Section IIIIBI Solid black lines with increasing 
thickness: i?g = 3000, 300, 30 EeV, and no,m = 0.8 nm.o.noovoi, 
and s = Snocvoi - 0.2, with nm.o.noevoi and Snoovoi as in- 
dicated in the second row of Table [TTl Black dashed hne: 
standard spectrum (Eq. I27p with uniform distribution of ini- 
tial voltages among magnetars, — 300 EeV, f2i,4 and 
?^m,o = 5.6 X 10~** Mpc^^ yr~^- AH cases are for a pure pro- 
ton injection, Si.min = 3 x 10^** eV, Si.max = 3 x 10^^ eV, 
fJ'33, and 771. The green dotted hne represents the 

LIGOIII and the Einstein Telescope (ET) approximate sen- 
sitivities [4^, [5^, the pink dotted lines the DECIGO and 
DECIGO Advanced sensitivities j6Q] and the purple dotted 
line the BBO sensitivity [sl, [5^, all in correlation modes. 
The values on the upper x axis represent the energy E = 



3 X 10^'Z77i/i33(vri^/10'' 



eV. 



[52l . |5J, |55| . The cross-correlation of two ground-based 
third generation LIGO-type interferometers (LlGOlIl or 
the Einstein Telescope) would lead to sensitivities inter- 
esting in the scope of this study [5^, ■ Sensitivities 
would be further increased with the subsequent genera- 
tion of experiments such as BBO and DECIGO, which 
should work in correlation mode by themselves, thanks 
to multiple detection units 58-60]. 

In what follows, we discuss the gravitational wave 
background spectra obtained for each distribution of 
magnetars discussed in Section IIIIl assuming that they 
are the sources of the observed UHECR flux. The ana- 
lytical calculations of these spectra are presented in Ap- 
pendix |D] The integrals are computed numerically. 

We present the cases of an SFR-type source evolution, 
which can be intuitively expected for a magnetar evolu- 
tion history, as magnetars are believed to form princi- 
pally in star-forming regions [s']. We conjecture that the 
parameters indicated in Table HTl in absence of evolution 
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FIG. 5: Red solid line: energy density of the stochastic grav- 
itational wave background figw produced by magnetars as a 
function of the observed frequency, assuming a distribution 
of /I among magnetars and a fixed v\ — fmax ~ 10* /vr, as in 
Section IIII CI We assume a pure proton injection, /? = 700, 
Si.min = 3 X IQi* eV, £i,n.ax = 3 X 10^1 eV, I45, and r,i. 
We take no^m = 0.8 nm.o.noovoi, and s = Snocvoi - 0.2, with 
Wm,o,noevoi aud Snoevoi as indicated in the third row of Ta- 
ble |lT] Black dashed line as in Fig. H) 



can be adapted for the SFR case by applying the rules 
stated at the end of Section IIIIDI 

In Figures 21 [S] and [51 the black dashed lines represent 
the gravitational wave spectrum obtained for the stan- 
dard magnetar population scenario, for which the initial 
voltages are uniformly distributed. It corresponds to the 
signal computed by Refs. 0, [5l|, though the normaliza- 
tion is different. As it appears in Eq. (P7)) . the quantity 
rigw scales as nm.Oj and it is hence straightforward to esti- 
mate the values that ilg^ could take for other source den- 
sities. In this study, we work under the assumption that 
magnetars are emitters of the observed UHECRs, which 
leads to the present scalings for gravitational waves. 

Figure 2] presents the gravitational wave spectra Jlgw 
obtained for different distributions of ili among magne- 
tars, with fj, fixed. The higher contribution of the low 
initial angular velocity sources (needed to fit the UHECR 
spectrum) enhances the gravitational wave spectrum at 
low frequencies, especially in the range i^o ~ 1 — 100 Hz 
that will be probed by future instruments in correlation 
mode. The poloidal field suggested from the observation 
of Galactic magnetars is of order ^ IQ-'^'' G, implying 
that e - IQ-^ - lO'^ (Sg > 3000) would be conservative 
values. The thin line of Fig. 2] shows that the gravita- 
tional wave signal would then still be a couple of orders 
of magnitude above the standard prediction, though be- 
low reach of the future instruments. The magnetar el- 
lipticity is however highly uncertain (see e.g., (6l| ) and 
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FIG. 6: Energy density of the stochastic gravitational wave 
background produced by magnetars as a function of the ob- 
served frequency, assuming Bd ~ avi (Eg. I16|) . with a — 10^^ 
and 10^'^ G s for the thin and thick lines respectively, and 
for [3 = 700 and 2000 for blue and red lines respectively. 
We assume a pure proton injection, -Ei.min = 3 x 10^* eV, 
iSi^max = 3 X 10^^ eV, 145, and 771. The parameters no,m 
and s are as indicated in the fourth row of Table In] with the 
adequate transformations to obtain the SFR evolution case. 
Black dashed line as in Fig. |4l 



more optimistic values of the critical gravitational en- 
ergy: _Eg = 30 and 300 EeV (corresponding to e ~ 0.3 
and 10~^ respectively) could be considered. The signal 
could then be strong enough to be detected by the future 
BBO and DECIGO, and could be close to detection for 
third generation LIGO-type instruments (LIGOIII or the 
Einstein Telescope). 

Note that the fact that a background appears to par- 
tially overlap to the sensitivity curve of a given detector 
does not necessarily imply that the signal is detectable. 
A calculation of the signal-to- noise ratio assuming a given 
observation time (as done in Refs. d^Ull) is needed to 
thoroughly evaluate the detectability. Such estimates are 
beyond the scope of this paper, in particular as the design 
of the instruments discussed here is still highly uncertain. 

In Figure [S] we plot the case where a distribution of 
/i is taken among magnetars, with a fixed initial angular 
velocity fij, and thus a fixed Vi = i^max- In this situa- 
tion, the dipole moment being independent of the grav- 
itational radiation frequency, we expect VL^^ to behave 
in the same way as for the uniform distribution case. 
The red dashed line in Fig. [5] arbors indeed the same 
shape as the one calculated in Eq. (P?]) . There is a factor 
~ (Aimin/Mmax)"'*Mmm(x9?77r2t'^a^/c2)(s - 1)"^ of differ- 
ence with the case of uniform distribution of voltages 
among magnetars (see Eq. ID6I) . This factor is smaller 
than unity for our choices of parameters. This case is 
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therefore "disapointing" in terms of gravitational wave 
signal, as it is expected to produce a signal weaker than 
in the standard magnetar population scenario. One may 
albeit point out that a sheer distribution of surface dipole 
field strengths _Bd, uncorrelated with the angular veloc- 
ities, might not be favored in the current magnetar for- 
mation models. 

Finally, Figure [5] shows an application for distributions 
of initial voltages among magnetars, with a dependence 
between the dipole magnetic field strength of the star 
and its initial angular velocity [/i = f{ili)]. We adopted 
Eq. as an example of such a relation, and computed 
the background gravitational wave spectrum for two val- 
ues of a. High values of a imply that high source den- 
sities are necessary to fit the UHECR spectrum and en- 
hance consequently the gravitational wave signal. For 
the most optimistic cases, the signal lies above BBO and 
DECIGO sensitivities, and is close to the sensitivity of 
third generation LIGO-type instruments. For low values 
of a however, the signal drops below the one emitted by a 
standard population of magnetars. The ellipticities rep- 
resented here are also optimistic (e ~ 10~^ for blue lines 
and 0.3 for red lines). More reasonable deformations of 
order e ~ 10"'^ would lead to a signal lower of about one 
order of magnitude. 



V. DISCUSSION AND CONCLUSION 

We have presented results for UHECR spectra pro- 
duced by newly born magnetars and their associated 
gravitational wave signatures assuming pure proton in- 
jection. The metal-rich supernova environment in which 
the magnetars form supports however the idea that a 
mixed or iron-rich cosmic-ray composition might be in- 
jected. As argued by Refs. [sl, [gI], strong electric fields 
could strip iron nuclei off the surface of the young neu- 
tron star. The magnetar magnetosphere could thus be 
enriched in iron-peak elements that can be accelerated 
via the unipolar induction mechanisms evoked in Sec- 
tion [Hi The escape of such accelerated nuclei is not ob- 
vious though: the rough estimates of Ref. [ll| show that 
ultrahigh-energy iron nuclei could survive the crossing 
of a supernova envelope under certain conditions, while 
Ref. [l2| proposes a different scenario in which it is sug- 
gested that only protons could escape. Given the recent 
results of Auger that could be interpreted as an indica- 
tion of a heavy composition at the highest energy end 
of the spectrum [6^, it will be of interest to investigate 
further the case of heavy nuclei. 

If a composition different from pure proton were 
injected, the spectral index necessary to fit the observed 
UHECR spectrum would be globally harder, implying 
that the distribution index of the potential drop among 
magnetars s should be changed accordingly (on average, 
one would need s' ^ s — 0.2 in each case). This 
discussion is valid also when different source evolutions 
are assumed: this should lead to a change in s, thus 



in a slight change in the normalization of the overall 
gravitational wave signal. The gravitational wave 
signature should then be slightly increased. On the 
other hand, the critical gravitational energy Eg scales 
as the particle charge Z: this should lead to a weaker 
gravitational wave production, all other parameters 
being fixed. Overall, the signal is not expected to be 
drastically modified in case of heavy nuclei injection. 

The single source injection spectrum calculated in 
Eq. [5] assumes a mono-energetic acceleration at a given 
time. Though this assumption is roughly valid for a pure 
electrostatic acceleration in the magnetar gap region, it 
might no longer stand for acceleration in the magnetar 
wind zone. More complex mechanisms than the toy 
induction model mentioned in Section should be at 
play. Wake-field acceleration or other stochastic accel- 
eration mechanisms as suggested by [H, l65l - [67j could 
lead to a natural spectrum. Other models such as 
magnetic reconnection models predict slopes in —1. The 
initial injection spectrum is at the moment uncertain 
and a thorough study of the acceleration mechanism 
that might be happening in the wind will be necessary. 
This paper demonstrates that even extreme injection 
cases can be reconciled with a proper distribution of 
parameters among the magnetar population. 

Figure [7| compares the gravitational wave backgrounds 
obtained in the present paper with the ones produced 
by other astrophysical or cosmological sources. A review 
on these backgrounds can be found in Ref. [s^. Unless 
the rates of binary neutron star coalescence in the 
Universe is over-estimated by many orders of magnitude, 
that signal should lie above the expected magnetar 
signals. The signal expected from r-mode instabilities 
in young fast-rotating neutron stars is less established, 
but under the optimistic assumptions taken in Fig. [71 
it could also lie above the magnetar spectrum. Most of 
the other astrophysical sources, including core-collapse 
supernovae, should produce weaker signals than the 
optimistic (though reasonable) magnetar populations 
invoked in this study. A good knowledge of the strength, 
frequency range, and statistical properties of overlapping 
backgrounds could help distinguish them from the 
signals that lie beneath. 

With improved statistics at the highest energies, and 
an increased anisotropy signal, one might be able to es- 
tablish if the sources of UHECRs are continuously emit- 
ting, or transient objects (see, e.g., [l3l)- In the latter 
case, long gamma-ray bursts (GRBs) and the birth of 
magnetars would be the best remaining candidates. The 
distinction between these two sources from their spatial 
distribution in the sky would be difficult: as transient 
sources, they should not be visible in the arrival direction 
of ultrahigh-energy events, and as they form in similar 
environments, they should present close spatial distribu- 
tions. 
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FIG. 7: Comparison with energy densities of the gravitational 
wave backgrounds produced by other astrophysical or cos- 
mological mechanisms. Black solid line: magnetar popula- 
tion with distribution of initial voltages assuming Bd = ai'i, 
with a = 10^^, red solid line: magnetar population with 
distribution of initial angular velocities f^i with /i fixed, for 
-Eg = 30 EeV, all the other parameters being as in Figs.|4]and 
|6l black dashed line: standard uniform case as in Fig. ID Back- 
grounds from binary neutron star coalescence (BNS) from 
[24| . r- modes assuming that 1% of newly born neutron stars 
cross the instability window from [s^], core-collapse super- 
nova from [6^. Horizontal dotted line: maximum version 
of the gravitational wave stochastic spectrum produced dur- 
ing slow-roll inflation assuming a ratio of the tensorial to the 
scalar contributions to the cosmic microwave background ra- 
diation anisotropy TjS = 0.3 and ±10~^ for the running of 
the tensorial power-law index [69l |. 



magnetars would no longer be an issue, but whether the 
acceleration happens by unipolar induction in the neu- 
tron star wind or by Fermi-type acceleration in the GRB 
shock waves would remain an open question, which could 
be investigated by the observation of gravitational waves. 
As we saw, the unipolar induction model requires the 
magnetar to have a certain level of pulsation and mag- 
netic fields that should naturally lead to the emission of 
gravitational waves. The detection of gravitational waves 
with the spectra that we predict here could be one evi- 
dence that high enough values for /i and Jli are reached 
and that unipolar induction acceleration should conse- 
quently take place. On the other hand, as mentioned 
before, the gamma-ray hurst in itself should only lead to 
a weak signal. 

Another secondary signature of UHECR acceleration 
in newly born magnetars that could also help distinguish 
magnetars from GRB scenarios is the neutrino emission. 
Ref. [gI] calculated that the neutrino flux produced by 
UHECRs crossing the magnetar surroundings could be 
detected by IceCube. This signal depends howbeit on 
the poorly known opacity inside the source (though 
the column density of the envelope the particles goes 
through is constrained by the fact that UHECRs need 
to escape, see Ref. [lH). The neutrino spectra found 
by Ref. [63 should also be affected by a change in the 
distribution of magnetar parameters to fit the observed 
UHECR spectrum. The spectra of cosmogenic neutrinos 
(produced by the interaction of UHECRs interacting 
with the radiative cosmic backgrounds while propa- 
gating from their source to the Earth) only depends 
on the overall UHECR spectrum shape, on the source 
evolution and the injected composition, see e.g., [36| and 
references therein. The distributions introduced here 
should thus not lead to new cosmogenic neutrino spectra. 



The detection of specific gravitational wave signa- 
tures might be one way of making this distinction. 
Long gamma-ray bursts are not believed to be strong 
gravitational- wave emitters, as the signal due to a burst 
(i.e., by collision of particles) should not be collimated 
and be too dilute to be observed by any instrument (e.g., 
[Tlj). It is possible that the progenitor of the GRB (e.g., 
collapsars, neutron star mergers, magnetar s...) pr oduces 
a certain level of gravitational radiation |72h74| |. The 
background spectrum created should broadly differ from 
the one predicted here for UHECRs. Moreover, the pa- 
rameters of the likely acceleration site in GRBs, the shock 
waves due to the explosion (t^, are not directly probed 
by gravitational waves from the central engine. The grav- 
itational radiation emitted by GRB progenitors might 
therefore not be connected to the capability of GRBs to 
accelerate UHECRs. 

There have been suggestions that magnetars could be 
the progenitors of long GRBs [76l - [78| . If this were 
the case, whether the sources of UHECRs are GRBs or 



We have demonstrated that, by relaxing the assump- 
tion of an uniform distribution of magnetar physical pa- 
rameters, it is possible to reconcile the spectrum of UHE- 
CRs produced by these objects with the observed data. 
The assumptions made on the magnetar population to 
best fit the observed UHECR spectrum leads to specific 
gravitational wave background signatures. In some mod- 
els for which the neutron star deformation by the inter- 
nal magnetic field is strong and the dipole surface field 
is a large fraction of the saturation field, the gravita- 
tional background signals overlap the sensitivity curves 
of satellites such as DECIGO or BBO, and are within 
reach of the sensitivity of the third generation LIGO- 
type detectors in correlation mode (LIGOIII, Einstein 
Telescope). The physical assumptions required for such 
a level of detection can be viewed as optimistic, but are 
still plausible in the magnetar formation scenarios. The 
signal obtained can be up to 3 — 4 orders of magnitude 
higher than standard predictions in the frequency range 
1 — 100 Hz that will be best measured by these future 
generation instruments. The detection of specific gravita- 
tional background spectra as we predict here could probe 
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newly born magnetars as UHECR accelerators, and help 
solve the long-standing question of the origin of these 
particles. 
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Appendix A: Cosmic ray spectrum produced by a population of magnetars 
with a distribution of initial angular velocities fli 

For a distribution of initial angular velocities ili with a fixed dipole moment fj, among the stars, one can write: 



dJ{E,EO _ 9 Ic' 



Integrating over Ei, one gets: 
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The hypergeometric function is noted: 
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Appendix B: Cosmic ray spectrum produced by a population of magnetars 
with a distribution of dipole moment /i 

We calculate the spectrum obtained for a distribution of dipole moments, when fij is fixed among the stars. Using 
Eq. ([15]), the critical gravitational energy can be expressed as a function of Ei as follows: 



2 /p;r<\l/2 



One can then write (using Eq. [3]) : 



dJ{E,Ei) 9 X p-i^ Ei \-^-\ Ei l + EEi/el \ 

- Wgcom-r^Icilin,^— E Tiass[E) — In — A ■ (^2) 



dEi "«'=°'"l67r^^^"''"'"i?,^ax '°'''^\Ei,n..J '''\El + Ef/el 
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The integration leads to (we note S = s + 1): 

9 

167r -E-ijnax L ^M.+ l-f^i It -C' > -C/i,min , 

where j^^^{E) and j^^+{E) are defined (we noted 5 = s + 1): 



Jm.-(^) 



-s 

i.max 



(1 - 1 



{l-S)ln[^i±^)+h{EE,/el)-2h'iE^/el)) 



(1 - S") In 



E el + El 
'Eiel + EEi 



E €^ + E: 



J J E 



The hypergeometric function h' is defined: 



, , 1 — s 3 — s 

h [x) = 2F1 1, — — , — — , -a; 



(B3) 

(B4) 

(B5) 
(B6) 

(B7) 



Appendix C: Cosmic ray spectrum produced by a population of magnetars 
with a distribution of dipole moment /i = f{Q,i) 

We assume a distribution of dipole moments with fi = /(J^i) defined in Eq. (jl6l) . For simplicity, we assume that Eg 
is fixed among magnetars (see discussion in Section fill D|) . Plugging in the relation obtained in Eq. ((T5)) between $i 
and ^ into Eq. (jAll) . one obtains: 



dJ{E, Ei) 



dEi 



9 Ic^ 



gcom „ T^uiX 1 J-^ 

167r q V^i.max 



A 



1/3 



Ei 



Ei, 



-s-1/3 



In 



i^i 1 + {E/Eg) 
E 1 + [E-jEg) 



where 



A — qrj- 



The integration is the same as in Appendix \K[ replacing the index by 5" = s + 1/3, which leads to: 



9 /c^ (A 

J{E) ^ W^goom-T rimX 

167r q V-E^i.niE 



1/3 



jn.-.S'{E) if £;<i;i, 
ja+,5'(^) if E>Ei, 



mm 
min 1 



where jfi^-^s'{E) and jn.+,S'{E) are the functions defined in Eqs. (jA4IA6P where we replace s by S' . 

Appendix D: Gravitational stochastic background spectra for various magnetar populations 

1. Distribution of Qi (thus of i/i) 



(CI) 



(C2) 



(C3) 



Assuming that the magnetic dipole moment does not vary from one source to another (i.e. remains constant), 
one can re-write the distribution of initial voltages (Eq. I12p as a function of the initial frequency Vi as follows: 



drini drim dEi 2qrjfi'K I 

avi dEi dui & 



Via 



(Dl) 



Taking into account the distribution of sources according to the initial voltage, i.e. to the initial frequency v\ under 
our hypothesis, Eq. (P7)) can be expressed 



f^gw(l^o) = 5.7 X 10 ( ) "m.OKl^Ql^i 



0.7 



2s 

i,max 



v\'''dv; I dz 



{i + zfn{z) di^ 



with 



2qrjinr'^ 



(D3) 
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2. Distribution of /i 

We assume that all magnetars have the same initial frequency Vi = t'max, and that jj, (thus i?d) varies from source 
to source. We can then write: 



In,,! dri„i d$i qr/Tr^ 2 / A* 



d^ d$i d^ 



(D4) 



And thus: 



/0.7 



17gw(i^o) = 5.7 X 10-'" ( — ) n,„,o 



■ 1^0 



dz — , ' ' , X 



(1 + Z)2f](z) 

-[1^0(1 + 



(D5) 



The analytical integration of the last integral over /i leads to: 



d/i 



gw 



dv 



- 1 



(D6) 



where h' is the hypergeometric function defined in Appendix [ 



3. Distribution of fj, — f{fli), i.e., Bd = ,f'(w) 

Under the assumption of Eq. (|16p connecting the dipole moment n to the initial angular velocity fii, the magnetar 
initial voltage can be expressed 

*. = ^ = ^-.'- m 

As a consequence, the distribution of magnetars according to the initial frequency reads: 



dz/i d$i S^-i d$i di/i 



(D8) 

nrnX^^ ^ ( ) . (D9) 

This yields the gravitational wave spectrum: 



f^gwM - 5.7 X 10-^6 ( I ^^ ^^-ll^lli::!^^ / ^M^) d^i X 



L 



,0 (1+Z)2f7(z) dv 

with 



^[(^0(1 + z), i^ij , (DIO) 



10 4o2nl0 2 



(D12) 
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